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1 Introduction and main result 

Consider an analytic system of the form 

x = y + X{x,y), y = Y{x,y), (1.1) 
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where X, Y = 0{\x, yp) for (x, y) near the origin. The following criterion for the existence 
of a center or a focus at the origin of (1.1) has been established in [U [T5] . 

Theorem 1.1. Let (1.1) have an isolated singular point at the origin. Let 

Y{x, F{x)) = ax'"-^ + 0(x2"), aj^O, 
BX BY 

— (x, F{x)) + — (x, F{x)) = + 0(x"), 

ax ay 

where y = F{x) is the solution of the equation y + X{x,y) = satisfying F{0) = 0. Then 
the origin of (1.1) is a center or a focus if and only if a is negative and + Aan < 0. 
Lyapunov [15j also introduced the generalized polar coordinates 

x = rCs{e), y = r''Sn{e) 

and the return map to give a way to find focal values in solving the center-focus problem 
for (1.1), where (Cs(t), Sn{t)) is the solution of the initial problem 

x = y, y = -x^^-\ (x(0),i/(0)) = (l,0). 

Sadovski [19] (see also [3J) and Moussu [16] investigated the problem using Lyapunov 
function (Lyapunov constants) and normal form, respectively. Then different ways of ob- 
taining the focal values, Lyapunov constants or their equivalent values and the bifurcation 
method of local limit cycles were further given by Chavarriga, Giacomini, Gine & Llibre 
[7], Alvarez & Gasull [L> 2J and Liu & Li [III [H [131 [S] • From Takens [21j we know that 
(1.1) can be formally transformed into a formal normal form 

x = y, y = -g{x) - yf{x), (1.2) 

where g{x) = ax"^ + 0{x'^~^^), m > 2 (the system (1.2) is a generalized Lienard system). 
Then, Strozyna & Zolfidek [20] proved that this formal normal form can be achieved 
through an analytic change of variables. Thus, if (1.1) has a center or focus at the origin, 
then it can be changed into (1.2) with 

g{x) = x2""i(a2„-i + 0(x)), n>2, a^n-i > 0. (1.3) 
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From Alvarez & GasuU [2] we see that under (1.3) through a change of variables x and t 
of the form 

r / X , , , , dt 

u = [2n / ^((xjaxj 2n (sgnx) = ti(a;), — 







dti g{x 
the system (1.2) can be changed into 



x = y, y = -x2"-i-y/», (1.4) 

where 



1.5) 



/» = x^^-^f{u-\x))/g{u-\x)), n > 2, 
u{x) = [2n g{x)dx]'^{sga.x) = (a2n-i)^(a; + 0(x^)). 
Then, by Theorem 1.1 system (1.4) has a center or a focus at the origin if and only if the 
function / given in (1.5) satisfies 

/»= J2 ^J-^'' d-4n<0. (1.6) 

j>n-l 

By Filippov's theorem (see e. g. Ye et al. [22]) under (1.6) the system (1.4) has a stable 
(unstable) focus at the origin if there exists an integer / with 21 > n — 1 such that 

b2l > Oi< 0), b2j = 0ioi J < I, (1.7) 

and it has a center at the origin if b2j = for all 2j > n — 1. 

Passing to the generalized polar coordinate {x, y) = {rCs{6),r'"'Sn{6)) we obtain from 
the system (1.4) the equation 



n+j 



The function on the right hand side of (1.8) is periodic of the period T = 2y^r(2^) /r(^^). 
Let r{6,ro) denote the solution of (1.8) with the initial value r(0) = tq. Then 

i>i 

Alvarez & GasuU [2] called the constant Vk the kth generalized Lyapunov constant of 
(1.8) assuming Vi = 1, V2 = ■ ■ ■ = Vk-i = 0. They also studied the normal form (1.4) and 
proved the following theorem. 

Theorem 1.2. Let (1.6) and (1.7) be satisfied. Then 



(2) Vi = 1, Vj = for 1 < j < 2 - n + 21, and V2_„+2/ = -Kih2i if either b^-i = 
or bn~-i 7^ and n is even, where Ki is a positive constant. 

For the case of n = 2 Liu and Li [11] introduced a different generalized polar coordi- 
nates of the form x = rcos6, y = r"^ sin 6 to change (1.1) into the form 

assuming the origin is a center or a focus. Let r{9, h) denote the solution of the 27r-periodic 

system 

dr R{e,r) 



de Qie,r) 

satisfying f(0) = h. Note that the initial value problem is well-defined also for negative 
h. 

For analytic functions 0, 0i, . . . , 0^ defined on a domain D we will write = O(|0i, . . . , 0^1) 
if there are analytic functions -01, . . . , -0^ on D, such that = '0i0i + ■ ■ ■ + ■0^0^ on D. 
Liu and Li [11] found the following facts. 

Theorem 1.3. Consider the system (1.1). Let the conditions of Theorem 1.1 be 
satisfied with n = 2 {or m = 3 ) such that the origin is a center or a focus. Then, 

(1) f{9,-r{n,h)) = -r{n-9,h); 

(2) A{h) = f(— 27r, h) — h = Yl,k>2 "^kh^ > where 

V2k+i = 0{\v2, ■ ■ ,V2k\), k>l; 

(3) the origin is a stable [unstable) focus if 

V2k < 0(> 0), and V2j = for j < k. 

In the latter case the origin is called a kth order weak focus of (1.1). 

Liu and Li ^llj also gave some new methods to compute the focus values V2, v^,--- , V2k, 
or equivalent values, and studied the problem of limit cycle bifurcations near the origin, 
finding a new phenomenon: a node can generate a limit cycle when its stability changes. 

In this paper we study the problem of limit cycle bifurcations near the origin for the 
analytic system 

x = y + X(x, y,5), y = Y{x, y, S), (1.9) 



where S = {6i, . . . ,6m) E D G M™" with D compact, and X, Y = 0(|x,?/p) for |x| small 
and 6 E D. Let y = F{x, 6) be the solution of the equation y + X{x, y, 6) = 0. We define 
the following two functions: 



^{x,F{x,6),6) + ^{x,F{x,6),6) 



g{x,6) = -Y{x,F{x,6),6), f{x,6) = 
By Theorem 1.1, if 

g{x,5) = x^^-\a2n-i{S) + 0{x)), n > 2, a2„-i(5) > 



(1.10) 



:i.iii 



fix,S)= J2b,{5)x\ bl_,iS)-Ana2n^iiS)<0, (1.12) 

j>n-l 

then the origin is a center or a focus of (1.9) for all S E D. 

Let us define a Poincare return map for the plane system (1.9). For each 6 E D 
and Xo 7^ with |xo| small consider the solution {x(t,xo,6),y(t,xo,6)) of (1.9) with the 
initial condition (x(0), y(0)) = (xq, F{xo, 6)). Then there is a unique least positive number 
T = t{xo, 6) > such that ?/(r, xq, 6) = F{x{t, xq, 6), 6) and Xox{t, xq, 6) > 0. See Figure 
1 for Xq > small. 















\ o 


A„ P{x„, 5)' 



Figure 1. The Poincare map of (1.9) with xq > 0. 

Thus, the Poincare return map is defined as 

, x{t,xo,5), < |xo| < eo, 

P{xo,6) 

"0, xo = 



where Sq is a small positive constant. Evidently, the function is uniquely defined, and 
it is continuous Bit Xq — under (1.11) and (1.12). Moreover, (1.9) has a periodic orbit 
near the origin if and only if the map has two fixed points near zero: one positive and the 
other one negative. For the analytical property of this function at xq = 0, we have the 
following theorem. 

Theorem 1.4. Let (1.9) satisfy (1.11) and (1.12) for all S E D. Then there is a 
unique analytic function P{xo, 6) in xq at xq = 0, satisfying ^(0, 5) > and having the 
expansion 



where P denotes the inverse of P in xq. 

Hence, the system (1.9) has a periodic orbit near the origin if and only if the analytic 
function d defined in (1.13) has two zeros in xq near xq = 0, among which one is positive 
and the other one is negative. The function d is called the succession function or the 
bifurcation function o/ (1.9) . 

The above theorem tells us that the function P{xo, 6) is analytic in xq at = as n 
is odd, and not analytic in xq at a^o = as n is even unless the origin is a center (in this 
case, P is the identity). 

For the property of the coefficients Vj in (1.13) we have further 

Theorem 1.5. Let (1.9) satisfy (1.11) and (1.12) for all 6 e D. Then 

(1) For n odd we have V2k — 0{\vi, v^, - ■ ■ , V2k-i\), k > 1. 

(2) For n even we have vi — 0, V2k+i — 0{\v2, v^, - ■ ■ , V2k\)-i k > 1. 

Define p„ = [1 + (— 1)"]/2. Then the conclusions of the above theorem can be written 
uniformly as 




(1.13) 



for \xq\ small, such that 



(1) if n is odd, then P{xo,S) — P{xo,5) for all \xo\ small; 

(2) if n is even, then for all \xo\ small 





for Xq > 0, 



V2k+p„ = 0{\vi+p^,V3+p„, ■ ■ ■ ,V2k-l+pJ), k>l. 
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From the proof of the above theorem we see that f 2A:+p„ depends on f , f 3+p„ , • • • , 
V2k-i+p„ smoothly. Using the theorem we derive the following two statements on limit 
cycle bifurcations near the origin. 

Theorem 1.6 (Bifurcation from Focus). Let (1.9) satisfy (1.11) and (1.12) for all 
5 E D. Denote p„ = [1 + (-l)"]/2. 

(1) // there is an integer k > 1 such that 

k+l 

l^2i-i+p„(5)| > for all 5 ED, 

i=i 

then there exists a neighborhood U of the origin such that (1.9) has at most k limit cycles 
in U for all 6 E D. 

(2) // there is 5q E D such that f2fc+i+p„(^o) 7^ 0, and 

V2j-i+pM = 0,j = !,■■■ ,k, 

then for an arbitrary sufficiently small neighborhood of the origin there are some 6 E D 
near 6o such that (1.9) has exactly k limit cycles in the neighborhood. 

Theorem 1.7 (Bifurcation from Center). Let (1.9) satisfy (1.11) and (1.12) for 
all 6 E D. Assume that there exist 6q E D and an integer k > 1 such that (1.14) is 
satisfied. If the origin is a center o/ (1.9) as V2j-i+p„{S) = 0,j = I,-- - ,k, then there 
exists a neighborhood U of the origin such that (1.9) has at most k — 1 limit cycles in 
U for all 6 E D near 6q, and also, for an arbitrary sufficiently small neighborhood of the 
origin there are some 6 E D near 6^ such that (1.9) has exactly k — 1 limit cycles in the 
neighborhood. 

The theorem means that the cyclicity of the system at the point 6o is equal to A; — 1. 
Now, different from [2] and [II]-[ll], we give the following new and more reasonable 
definition. 

Definition 1.1. We call f2fc+i+p„(5) the generalized focal values of order k of (1.9) 
at the origin. 

By Theorem 1.6, we see that a nilpotent focus of order k generates at most k limit 
cycles under perturbations as long as the perturbations always satisfy (1.11) and (1.12). 
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The generalized focal values Ui+p^ , Us+p^ , • • • ,V2k+i+p„, ■ ' ' can be calculated using 
the normal form of system (1.9). We will give a method how to do it. By Strozyna and 
Zol§-dek [20] we know that (1.9) has the following analytic normal form: 

x = y, y= -9{x, 6) - yf{x, 6). (1.15) 

We remark that here / and g in (1.15) may be different from ones given by (1.10). As 
before, let 5 G -D C M™ with D compact. Also, suppose for small the function g{x,6) 
satisfies (1.11). Define 



F{x,6) = / f{x,6)dx, G{x,6) = / g{x,6)dx. 
Jo Jo 

It is easy to see that the equation G{x, 6) = G{y, 6) for xy < defines a unique analytic 

function y = a{x, 6) = —x + O(x^). Introduce 

F{a{x,6),6) - F{x,6) = ^Bj{6)x^. (1.16) 

i>i 

By Theorem 1.1, if (1.15) satisfies (1.11) and (1.12) then it has a center or focus 
at the origin. Thus, under (1.11) and (1.12) the Poincare return map for (1.15) is well 
defined near the origin. 

Theorem 1.8. Let (1.15) satisfy (1.11) and (1.12) for all 5 e D. Then for xq > 
small, the Poincare return map P{xq, 6) has the form 

j>0 

where P*{xo,6) = 0{xo), 

^i+p„(5) = K*B2l+^i6) + (1 -p„)0(i?2Vi), 

V2j+l+p„{6) = Kl_^jB2l+2j+l{S) + 0{\B2l+l, B2I+3, • • • , B2l+2j-l\), j > 1, 

/ = [n/2], and K^^j, j > are positive constants. Thus, Theorems 1.6 and 1.7 hold if 
V2j+i+p„ IS replaced by B2i+2j+i, J > 0. 
Let 

f{x,6) = J2bM^'- (1-18) 

i>o 

Then, we further have for (1.15) 



Theorem 1.9. Let (1.15) satisfy (1.11), (1.16) and (1.18) for all S E D. Assume 
there exist 5q E D and k > [n/2] such that 

52fe+i(<5o) <0(>0), 52,_i(5o)=0,j = l,---,A;. (1.19) 

Let one of the following conditions be satisfied: 

(a) n = 2, and 

boiSo) = 0, bliSo) - Sa^iSo) < 0; (1.20) 

(b) n>2, g{-x,6) = -gix,5), f{-x,5) = f{x,6), and 

b.iSo) = for J = 0,--- ,n-2 and d(5o) - 4na2„^i(5o) < 0. (1.21) 
Then we have 

(1) For 6 = 6o (1-15) has a stable (unstable) focus at the origin. 

(2) If further 

d{Bi, B-i, - ■ ■ , B2k-i) ,^ s , 

rank——— ——[do) = /c, 

d{6i,62,--- ,Sm) 



then for an arbitrary sufficiently small neighborhood of the origin there are some S E D 
near 5 such that (1.15) has at least k limit cycles in the neighborhood. 

From Theorems 1.4-1.8, it seems that under (1.11) and (1.12) we have solved the 
problem of limit cycle bifurcation for generic systems. Theoretically it is, but in practice 
it is not. The reason is that in general we do not know what is the transformation from 
(1.9) to its normal form (1.15). Here we give a method to solve the problem completely 
both theoretically and in practice. It includes three main steps below. 

First, under (1.11) and (1.12) by the normal form theory (see, for instance, [21]), for 
any integer m > 2n — 1 there is a change of variables of the form 




+ Hmiu,v,5), 



X 

y 

where Hm{u,v,6) = 0{\u,v\) is a polynomial in u,v of degree at most m, such that it 
transforms (1.9) into (1.22) (called the normal form of order m of (1.9), or the Takens 
normal form; we still use [x, y) for the new variables u, v) 

x = y + Xm+i{x, y,5), y = -g^^x, 5) - yfm-~i{x, 5) + Y^+i{x, y, 5), (1.22) 



where 

m m—1 
j=2n—l j=n—l 

with a2n-i{S) > and bl_i{S) - Ana2n-i{S) < 0, and X^+i{x,y,S), Y^+i{x,y,S) being 
analytic functions satisfying X^^i, l^+i = 0{\x,y\'^'^^). Here, we should mention that 
the functions Qm and fm-i depend only on the terms of degree at most m of the expansions 
of the functions X and Y in (1.9) at the origin. 

The Poincare maps of (1.9) and (1.22) are essentially the same. We can suppose that 
the Poincare map of (1-22) is P{xq, 6) having the expansion 

P(xo, S)-xo = J]^^i(5)4 (1-23) 
j>i 

for xq > small. 

Second, truncating the higher order terms in (1.22) we obtain the following polyno- 
mial system of degree m 

i = y, y = -gm{x, 5) - yfm-i{x, 5). (1.24) 

In practice, for given system (1.9) it is not difficult to find the corresponding system 
(1.24). For (1.24) we can further use Theorem 1.8 to find its focal values at the origin up 
to any large order. Let Pm{xo: S) denote the Poincare map of (1.24). It has the expansion 

Pm{xO:S) - Xo = ^Vj{5)xi (1.25) 

for xo > small. 

Third, we want to use Vj{6) for Vj{6). Here, a problem we would like to solve is the 
following: For any given k > 1 find m > 2n — 1 such that Vj{S) — Vj{S) for 1 < j < k. 
The following theorem gives an answer. 

Theorem 1.10. Consider (1.22) and (1.24). Then for any integer k > 1, if m > 
{k + 2)n - 2 then 

Vj{S) = Vj{S) for l<j<kn. (1.26) 

Therefore, we have 
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Corollary 1.1. Under (1.11) and (1.12) for any integer k > 1 for (1.9) the coeffi- 
cients f 1, f2, ■ ■ ■ , ffcn in (1.13) depend only on the terms of degree at most {k + 2)n — 2 of 
the expansions of the functions X and Y at the origin. 

Obviously, in the case of n = 1 (the elementary case), the above conclusion is a 
well-known results. 

We organize the paper as follows. In section 2 we first give preliminary lemmas. In 
section 3 we prove our main results. In section 4 we provide some application examples. 



2 Preliminaries 

Consider (1.9). In this section we will always suppose that (1.11) and (1.12) are satisfied. 
Introducing a new variable v = y — F{x, 6) we can obtain from (1.9) (reusing y for v) 

X = y{l + Zi{x,y,6)), 

(2.1) 

y = -9{x, S) - yf{x, 6) + y'^Z2{x, y, S), 

where the functions / and g are given by (1.10), and Zi and Z2 are analytic functions 
near the origin with Zi{x,y,6) = 0{\x,y\). In the discussion below we will often omit 6 
for convenience. As in Liu and Li [H] we will make a change of variables to (2.1) using 
the generalized polar coordinates 

x = rcosd, y = r"'sin6, r > 0. (2.2) 

Lemma 2.1. Let (1.11) and (1.12) be satisfied. Then the transformation (2.2) 
carries (2.1) into the form 



e = s{e,r) = ^^[So{e) + o{r)], 

r = R{e,r) = j^^[Roie) + Oir)], 



(2.3) 



where S and R are 2n-periodic in 9, and satisfy 

Sin + (-1)"-^^^, -r) = (-l)"-^^(^, r), R{n + (-1)""^^^, -r) = -i?(^, r), (2.4) 

and H{0) = cos^ 6 + n sin^ 6 > 0, 

So{e) = -[n sin^ O + bn-i cos" 6 sin 6 + a2n-i cos^" 0] < 0, 
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i?oW = COS ^ sin ^(1 - a2n-i cos^""^ ^ - 6„-isin^cos"-2^). 



Proof. From (2.2) we have 

X — cos Of — r sin y = nr""-*^ sin Of + r" cos 00. 

We solve the above equations for and r, and obtain (2.3) with 

cos Oy — nr^~^ sin Ox 



S{0,r) 



R{0,r) 



r"(cos^O + nsm'^0) ' 
sin Oil + r""^ cos Ox 



r"-i(cos2^ + nsin^ 0)' 
Then noting that 

cos(7r ±0) = — cos 6*, sin(7r ± 6*) = =F sin 6* 

and that (2.2) is invariant as {0, r) is replaced by (tt + (—1)""^^, — r) one can prove (2.4) 
easily. The other conclusions are direct. This ends the proof. 

By (2.3) and (2.4) we obtain the following analytic 27r-periodic equation 

^=mr), (2.5) 

where 

. sin ^^?/ + r"^^ cos 6'i; 

R{0,r) =r 



cos 9y — nr"- ^ sin Ox 
= r[Ro{0)/S,{e)+0{r% (2-6) 

^(tt + (-1)"-^^, -r) = {-1YR{0, r). 

Let r{0,h) denote the solution of (2.5) with the initial value r{Qi) — h. For properties 
of the solution we have 

Lemma 2.2. The solution r{0,h) = 0{h) is analytic in {0,h) for \h\ small, and 
satisfies 

(1) r(0, -r(n, h)) = -r(n + (-1)"-^^, h); 

(2) r{9±2n,h) = r{e,r{±2n, h)) . 

Proof. Let f{0) = -r(7r + (-l)'*"^^, h). Then by (2.5) and (2.6) we have 

f = (-l)"i?(7r + (-1)"-^^, r{n + {-1)^-^0, h)) 
= (-l)"i?(7r + (-l)*^-^^, -f{0)) 
^R{0,f{0)). 
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This means that f{9) is also a solution of (2.5). Then the first conclusion follows by the 
uniqueness of initial problem. The second one follows in the same way. This completes 
the proof. 

Further we have 

Lemma 2.3. Let P{xo,S) be the Poincare return map of (1.9) defined in section 
1. Then for \xq\ > small we have P{xo,5) — r(— 27r,xo) for xq > 0, and P{xq,5) — 
r((-l)"27r, xo) for xo < 0. 

Proof. First, it is easy to see that (1.9) and (2.1) have the same Poincare return 
map P{xo,S). Then, noting that ^ < for r > small by (2.3), by the definition of P 
and (2.2) we can see that 

P{xo, S) = x{t, Xo) = r(-27r, xo) 

for xq > small. Now consider the case of xq < 0. Let r*{9,h) denote the solution of 
(2.5) satisfying r*(7r) = h. Then we have similarly 

P{xo,5) =x{t,xo) = -r*{-n,-xo), 

since under (2.2) the points (xq, 0) and {P{xo, S), 0) on the {x, |/)-plane correspond to the 
points (tt, —Xq) and (— tt, —P{xo,6)) on the (6',r)-plane respectively. 
Further, by Lemma 2.2(1) we have 

r*{6, —h) — —r{TT — 6, h) for n even, (2.7) 

and 

r*{e, -r(27r, h)) = -r(7r + 9, h) for n odd. (2.8) 

Noting that by Lemma 2.2(2), xq = r{2TT,h) if and only if /i = r(— 27r, a;o), we see that 
(2.8) becomes 

r*(e, -Xo) = -r(7r + 6, r(-27r, xo)) for n odd. (2.9) 
Therefore, for < by (2.7) and (2.9) 

{r(2Tc,xo) for n even, 
r{—2n,xo) for n odd. 
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This ends the proof. 

Lemma 2.4. Let d{xo,5) — P{xo,5) — Xq. Then there exists an analytic function 
K{h, S) for \h\ small with K{0, S) = ^(tt, 0) > such that 

d{xo, S) = -K{xo, S)d{xo, d) (2.10) 

for Xq > small, where Xq — — r(7r, Xq). 
Proof. By Lemma 2.2, we have 

r((-l)"27r,xo) = -r{-7i,xo) = -r(7r, r(-27r, xq)). 
Hence, by Lemma 2.3 for xq > 

d(xo,5) = r((-l)''27r,xo) - xo 

= -r(7r, r(-27r, xq)) + r{n, xq) 
= -K{xo,S)[r{-27r,xo) - xq] 
= -K{xo,5)d{xo,S), 

where 

dr 

K{xo,S)^ / - — {n, Xq + s{r {-27:, Xq) - xo))ds. 

Jo 0*^0 

It is obvious that K is analytic for |a;o| small and K{0, 6) = ^(tt, 0) > 0. This completes 



the proof. 



3 Proof of the main results 

In this section we prove our main results presented in Theorems 1.4-1.10. 

Proof of Theorem 1.4. We take P{xo,S) — r(— 27r, Xq) for \xo\ small. Then by 
Lemma 2.2 P is analytic. Note that by Lemma 2.2, r(27r, xq) is the inverse of r(— 27r,xo) 
in Xq. Then Theorem 1.4 follows directly from Lemma 2.3. The proof is complete. 

Proof of Theorem 1.5. There are two cases to consider separately. 

Case A: n odd. By (1.13) and Theorem 1.4(1), we have 

d{xo, S) = d{xo, S) = ^Vj{S)x^o (3.1) 
14 



for all l^ol small. 

By Lemma 2.4, we can suppose 

K{xo, 5) = ^ kjX^, xq = -r(7r, xq) = ^ IjX^, (3.2) 

j>0 j>l 

where ko > 0, li — —ko. Substituting (3.1) and (3.2) into (2.10), we obtain 

j>l i>l i>0,j>l 

Comparing the coefficients of the terms Xq, Xq and x'fj' on both sides yields 

V2II + Vil2^ -(koV2 + kiVi), 

vjj + 3V3III2 + ^2(^2 + 2/1/3) + ViU = -{koV4 + kiVs + k2V2 + ksVi), 



V2j 



+ l'2j-livlj(Zl, ^2) + • • • + V2L2j-2,j{h, h,--- , hj-l) + Vil2j — — Yl^Lo^ kiV2j-i, 



where Lij{li, h, - ' ' ? ^i+i)) ^ = 1? 2, • • • , 2j' — 2, are all polynomials. Then from the above 
equations we obtain 

V2k = 0{\Vi,Vs,- ■ ■ ,V2k-l\), k>l. 

Case B: n even. By (1.13) and Xo — P"^(P(xo, S), S) we can find 

P~^{xo, S) = vixo + V2xl + V3XI H , (3.3) 



where 



^1 = {vi + 

V2 = -V2{V1 + 1)~^, 



(3.4) 



= -'"ji'^l + 1) ^^^^^ + Lj{v2, Vs,--- , Vj-i), 



where each Lj is a polynomial of degree at least 2. Now we suppose Xq > 0. Then (3.1) 
holds by Theorem 1.4. Further, noting that 5;o < by Theorem 1.4 again 

d{xo, S) = P{xq, 5)-Xo = P~^{xo, S) - Xq. (3.5) 
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Then, inserting (3.1), (3.2), (3.3) and (3.5) into (2.10) we obtain 

{vi - l)/i = -kovi, 

{vi - l)l2 + V2II = -{koV2 + kiVi), 

(3.6) 

{vi - + Lj{v2, V3, - ■■ , Vj-i) + Vjl{ = -{koVj + kiVj-i H h fcj-iVi), 

where 

Lj{v2,V3, ■ ■ ■ ,Vj-i) e (^2,^3, • • • ,Vj-l)- 

Finally, noting that li — —ko and substituting (3.4) into (3.6) we easily see that 

V2j+l = 0(\V2, V4,-- - , V2j\), j > 1. 

This ends the proof. 

Proof of Theorem 1.6. For the first part, suppose the conclusion is not true. 
Then there exists a sequence {dm} in D such that for 5 = 5^ (1.9) has k + 1 limit cycles 
Lm,i, Lm,2, ■ ■ ■ , Lm,k+i which approach the origin as m — > 00. Then by Theorem 1.5, the 
function d{xo, Sm) has 2A; + 2 non-zero roots in xo which approach zero as m — > 00. 

Since D is compact, we can assume 6m E D a.s m ^ 00. By our assumption, 

E,^;; \v2j-i+pM\ > 0. Thus, for some l<l<k + l, 

V2i-i+p„{So) 7^ 0, V2j-i+p^{So) = for 1 < j < / - 1. 

Therefore, by (1.13) and Theorem 1.5, we have 

d{xo, So) = V2i-i+p„{do)xl'-'^^" + 0{xl'^^"). 

Note that d{0, 6) = 0. It follows from RoUe's theorem that for some > the function 
d{xo, S) has at most 21 — 2+pn non-zero roots in (—£0, £0) for all |(5 — (5o| < £o- We have 
proved that the function d{xo, 5m) has 2A;-|-2 non-zero roots which approach zero as m — > 
00. It then follows that 2A; + 2 < 2l — 2+pn, contradicting to 2/ — 2+p„ < 2k+Pn < 2k + l. 
The first conclusion follows. 
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For the second one, by Theorem 1.5, the function d can be written as 

rf(xo,5) = 5^t;2,_i+p„(5)^o'"'^'"(l + ^.(a;o,5)), (3.7) 
i>i 

where Pj(0,5) = 0. Like in [8] one can show that Pj are series convergent in a neighbor- 
hood of 5q (see also e.g. [18l[T7]). Further, by (1.14), we can take f i+p^, f3+p„, ■ ■ ■ , f2fc-i+p„ 
as free parameters, varying near zero. Precisely, if we change them such that 

< l^l+Pnl < l^3+p„| < ■ ■ ■ < \v2k-l+p,, \ < 1, Vi+p^V^+p^ < 0, ■ ■ ■ ,V2k-l+p„V2k+l+p„ < 0, 

then by (3.7) the function d has exactly k positive zeros in xq near xq = 0, which give k 
limit cycles. This finishes the proof. 

By Theorem 1.4 and (3.7) we immediately have 

Corollary 3.1. Let (1.9) satisfy (1.11) and (1.12) for a fixed 6 eD. Then, if 
V2k+i+pS^) < 0(> 0), V2j-i+pS^) = for j = 1, ■ ■ ■ , /c 
the origin is a stable [unstable) focus of order k of (1.9). // 

f2i-i+p„(^) = for all j > 1 

the origin is a center of (1.9). 

Proof of Theorem 1.7. Under (1.14) t>i+p^ , f 3_|_p„ , ■ ■ ■ , V2k-i+p„ can be taken as free 
parameters. Further, by our assumption, the origin is a center of (1.9) as V2j~i+p„{S) = 
0, j = 1, ■ ■ ■ , k. It then follows 

V2j-i+p„{S) = 0(|wi+p„,U3+p„, ■ • • ,W2fc_i+pJ) for all j >k + l. 

Therefore, (3.7) can be further written in the form 

k 

d{xo,6) = 5^t;2,„i+p„(5)x^^"^+^"(l + P,(a;o,5)), 
i=i 

where Pj{0,5) = and Pj are series convergent in a neighborhood of 5q ([5]). Using the 
reasoning of Bautin [H] (see also e.g. [HI [13, [IB]) one can easily see that the conclusion of 
the theorem holds. The proof is completed. 
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Proof of Theorem 1.8. Now we consider (1.15), where g satisfies (l-H). Let 

F[x,5) = i f{x,5)dx, G{x,5) = I g{x,5)dx. 
Jo Jo 

If / satisfies (1-12), then the origin is a center or focus of (1.15), and 

F(a{x, 5), 5)- F(x, S) = J2 Bj(S)x^ = J2 (^.S) 

j>n j>ni 



where 

f— 1~)" — 1 rn" 

Bn^- 6n-i, ni^2l + l,l= - 

n 121 

and a{x, 5) — —x -\- 0{x^) satisfies G{a{x, S),6) — G{x, 5) for j^l small. Note that (1.15) 
is equivalent to the following system 

x^y- F{x, S), -g{x, 5) (3.9) 

which has the same Poincare return map P{xo,S) as (1.15). Introducing the change of 
variables x and t 



[2nG{x,S)]^n[sgnx) = (asn-i)^" (x + 0(a;^)) = (p{x), 



dti g{x,S) 
the system (3.9) becomes 

u^y-F{u,S), y^-u'^-\ (3.10) 

which is equivalent to 

u = y, y = -u'--^-yf{u,d), (3.11) 

where 

dF 

F{u,S) = F{^-\u),S), f{u,6) = ^(«,5)- 

The systems (3.10) and (3.11) have the same Poincare return map, denoted by Pi{uo, 5). 
One can see that the maps P and Pi have the relation Pio (p — (po P. Hence, 

P{xo, 5) -xo^ K{uo){Pi{uo, 5) - uq), 

where K{uo) — (a2n-i)~^ + 0{uo) is analytic. By (1.13) and (3.7), for uq > small we 
have 

Pi{uo,5)-uo = Y.'^ASX = $^t;2,-i+p„(5)«^^-'+^"(l + P,(uo,5))- 
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Hence, 

P{Xo, 5)-Xo = J2j>l V2j-l+pA^){a2n-l)~^Uo'^^''"il + Pjiuo, 6)) 

where Pj{uo,6) = 0{uo), P*{xo,6) = 0{xo). 

Since a satisfies G{a{x,S),S) — G{x,d) and xa < for \x\ small, we have (f{a) — 
—<f{x) or a — (fi''^{—(fi{x)) — (fi~^{—u), where u — (fi{x). Thus, we have 

F{a{x, 5), 5)- F(x, 5) = F(ip-\-u), 5) - F(ip-\u), 5) = F(-u, 5) - F{u, 5). (3.13) 

Let 

f{u,s)= mu^- 

j>n-l 



Then 



Thus, by (3.13) we have 



j>n ^ 



F{a{x,5),5)-F{x,S)^-2 ^ 



j>[n/2] 

Substituting u = (p{x) = (a2n-i)^(a; + 0{x'^)) into the equality above and comparing 
with (3.8) we obtain 

B21+1 = -Kihi, B21+2 = 0(b2i), 

B2l+2j+l = -Ki+jb2l+2j + 0{\b2U b2l+2: " " ' , 62Z+2j-2|), (3-14) 
-B2i+2j+2 = 0{\b2l,b2l+2, ■ ■ ■ ,b2l+2j\), j > 1, 

where Ki, Kij^i, ■ ■ ■ are positive constants. 

Then by Theorem 1.4 for i^o > small we clearly have 

Pi{uo, 5) = -uo + ^ Vj{5)u^Q = ^ yA^Wo, 

where Vj are introduced before Theorem 1.2. Thus, by Theorem 1.2, we have 

^"1 = -Kib2i + 0(bli), V2j+i\vi=-=v2j-i=o = -Ki+jb2i+2j, J > 1 
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(3.15) 



(3.16) 



for n = 2/ + 1 odd, and 

V2 = -Kib2U V2j+2\v2 = ---=V2j=0 = --^/+j^2«+2i , j > 1 

for n = 21 even, where -ft'i+j, J > are positive constants. Hence, 
wi+Pn = -Kihi + (1 -p„)0(6|), 

^2j+l+p„ = -Ki+jb2l+2j + <^{b2l,b2l+2, ■ ■ ■ ,b2l+2j-2), j > 1, 

where <^(0, 0, ■ ■ ■ , 0) = 0. Note that (1.15) is analytic in each bj. It follows from Theorem 
1.4 that is analytic in (62;, &2i+2, ■ ■ ■ , &2«+2i-2), which yields = O^hi, 62/+2, ■ ■ ■ , &2«+2j-2|)- 
Then (3.14) and (3.15) together give 

^l+p„ = §52m + (l-Pn)0(i?2VJ, 

V2j+l+p„ = J^B2l+2j+l + 0{\B2l+l, B2I+3, ■ ■ ■ , B2l+2j-l\), j > 1, 

Then (1.17) follows from (3.12) and (3.16). This finishes the proof. 

Proof of Theorem 1.9. Let |(5 — 5o| be small. For n = 2 we have j9„ = 1. Then 
the first conclusion is direct from Corollary 3.1 and (1.17)-(1.20). In fact, we have by 
Theorem 1.8 

V2ki5o) = KlB2k+iiSo), Kl > 0, V2j{5o) = for j = 1, ■ ■ ■ , A; - 1. 

For the second conclusion, we first keep Bi{6) = 0, and vary B-^IS), - ■ ■ , -62^-1 (^) near 
zero to obtain exactly k — 1 simple limit cycles near the origin. These limit cycles are 
bifurcated by changing the stability of the focus at the origin k — 1 times. Then we vary 
Bi such that < ^ |i?3|, and -B1-B3 < 0. This step produces one more limit cycle 
bifurcated from the origin by changing the stability of the origin which is a node now by 
|10j . The theorem is proved for the case of n = 2. 

For n > 2 since g{—x, 6) = —g{x, 6), f{—x, 6) = f{x, 6) we have 

bj{6) = for j = 0, ■ ■ ■ ,n-2 and &^_i(5) - 4na2n-i(5) < 

if 

52;+i(5) < 0(> 0), 52,_i(5) = 0,j = l,--- ,Z (3.17) 
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for some [n/2\ < I < k. In this case the origin is a stable (unstable) focus of (3.9) by 
Theorem 1.8. If (3.17) holds for some < / < [n/2\, then by [lOj again the origin is a 
stable (unstable) node of (3.9). Then the proof in this case is just similar to the above. 
This finishes the proof. 

We remark that if g{—x, S) = —g{x, 6) then a{x, S) = —x. 

Proof of Theorem 1.10. Consider (1.22). Without loss of generality, we can 
assume X^+i = in (1.22). Otherwise, it needs only to introduce a change of variables 
V = y + Xm+i{x,y). In this case, we can write (1.22) into the form 

x = y, y = -g{x) - f{x)y + y^^ipj{x)y\ (3.18) 

where 

g{x) = gUx) + 0{\xr+'), fix) = /^_i(x) + 0{\xn, = Odxp-^"^). (3.19) 

For the sake of convenience below, we rewrite the functions g, f and (pj as follows: 

g{x) = a;^"^-'^[(7o(x) + x'^gi{x) + x'^'^g2{x) + ■■■], 

/(x) =x"-i[/o(x)+x"/i(x)+a;2"/2(x) + ■■■], (3.20) 

ipj{x) = (fjoix) + x''-'^ipji{x) + x2""Vi2(a;) H , 

where fj, gj and ipji, j > 0, / > 1, are polynomials in x of degree at most n — 1, and ipjo, 
j > 0, are polynomials in x with degree at most n — 2. 

Now we change (3.18) by using (2.2) to obtain (2.5) satisfying (2.6) where 

X = r" sin^, y = r^""^ ^ Vj{e, r)r^'^, 

and by (3.20) 

Vo{9, r) = — cos^"""*^ 6'(7o(r cos 9) — sin" Ofoi^r cos 9) + r sin^ 9(poQ{r cos 9), 

Vj{9,r) = -cos2"-i+-'" %(rcos^) - sin" ^ cos-'" ^/j (r cos^) + r sin^^^' %o(r cos ^) 
i-i 

+ J2 sin^^'' cos(^'~'^)"^i 9^k,j-k{r cos 9), j> 1. 

A:=0 

(3.21) 

Hence, we obtain from (2.5) 
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where 

So{9,r) = -nsiry^e + cose Vo{e,r), Ro{e,r) = sin ^ cos ^ + sin ^ Vo{9,r), 

^3.22^ 

Sj{9,r) = cos 9 Vj{9,r), Rj{9,r) = sin 9 Vj{9,r), j > 1. 
By (3.21) and (3.22) we can further expand Sj and Rj in r to obtain for j >0 

n—l n—1 

SA9,r) = Y,Si+jn{9y, R,{9,r) = Y,Ri+j^{9y (3.23), 

1=0 1=0 

so that the above differential equation can be written as 



Further, letting 



and 



d9 Zj>oSj{9y 



we obtain 



Rj{9) = J2 M0)Si{9), J > (3.24) 

k+l=j 

^-rJ^RA&y- (3-25) 

j>o 

Note that for any j > 0, Sj depends only on with < k < j. Then by (3.24) one can 
see that 

For any j > 0, Rj depends only on R^ and Sk with < k < j. (3.26) 

Let r{9,ro) denote the solution of (3.25) with the initial value tq. The for ro small 
we have 

^(^^'^o) = ^rj{9yo 

where ri, r2, ^3, • • • satisfy ri(0) = 1, r2(0) = r3(0) = • • • = 0, and 

r[ = Ron, 

r'2 = Ror2 + Rirj, 

r's = R-oTs + 2Ririr2 + -^2^1, 
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which imphes that for any j > 1, the function Vj depends only on with < k < j — 1. 
Hence, by Lemma 2.3, (1.23) and (3.26) we come to the following conclusion: 

For any j > 1, Vj{d) depends only on and Sk with < k < j — 1. (3.27) 

Further, by (3.21)-(3.23), one can observe that for < / < n — 1, and Ri depend 
only on the coefficients of degree / of the polynomials go, fo and xipoo in x. Hence, by 
(3.27) we see that for 1 < j < n, Vj depends only on the coefficients of degree at most 
J — 1 of the polynomials go, fo and xtpoo in x. 

Similarly, for j > I and < / < n — 1 or jn < / + jn < (j + l)n — 1, Si+jn and Ri+jn 
depend only on the coefficients of degree / of the polynomials gj, fj, x(fijo and (fiij-i with 
i — 0, ■ ■ • , j — 1 in X. In other words, for jn + 1 < u < {j + l)n, Su-i and Ru-i depend 
only on the coefficients of degree u — 1 — jn of the polynomials gj, fj, XLpjo and ^ij-i 
with i = 0, • • • , j — 1 in a;. Let N]^a,h] denote the set of integers in the interval [a, b]. Then 
for jn -\- 1 < u < {j -\- l)n, we have 

i-i 

N[o,u-l\ = U ^[m,(i+l)n-l] (J-/V[jn,u-l]. 

i=0 

Thus, for all k e iV[j„^(i+i)„_i], Sk and Rk depend only on gi, fi, xipio and Lpi-i-i with 
I — 0, • • • ,i — 1. And for k e N^n^^^i^, Sk and Rk depend only on the coefficients of 
degree k — jn of the polynomials gj, fj, x(fjo and (fii,j-i with I — 0, ■ ■ ■ , j — 1 in x. 

Therefore, by (3.27) for jn + 1 < u < {j + l)n, f„(5) depends only on the functions 
gi, fi, xipio and (pi^i-i with / = 0, • • • ,i — 1, i = 0, ■ ■ ■ ,j — I and the coefficients of degree 
at most u — 1 — jn of the polynomials gj, fj, x(fijo and fij-i with I — 0, - ■ ■ , j — 1 in x. 

We claim that if j > 0, m > {j + l)n, then for jn + 1 < u < {j + l)n, Vu{S) depends 
only on the functions gi, fi, with i = 0, - ■ ■ ,j — I and the coefficients of degree at most 
u — 1 — jn of the polynomials gj, fj in x. 

In fact, by the above discussion, we need only to prove (foo = in the case j — and 
ipis — for I + s < j and < / < j — 1 in the case j > 0. This can be shown easily since 

ipjo = 0(|xp-^-^), ifj, = 0(|a;p-^'-"") for s > 1 

and 

deg (fijo < n — 2, deg (fjs < n — 1 for s > 1 
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by (3.19) and (3.20). 

By (3.20) again, the above claim can be restated that if j > 0, m > (j + l)n, then for 
jn+1 < u < {j + l)n, Vu{S) depends only on the coefficients of degree at most 2n + u — 2 
of g and the coefficients of degree at most n + -u — 2 of / in x. Thus, for any integers k 
and m satisfying k > 1 and m > {k + l)n, by taking j = 0, ■ ■ ■ , k we know that for all 
1 < u < {k + l)n, Vu{S) depends only on the coefficients of degree at most 2?7, + m — 2 of 
g and the coefficients of degree at most n + u — 2 of f in x. 

Finally, by (3.19), if m > {k + 3)n - 2 then 

2n + u ~ 2 < m, n + u — 2<m — 1 for u < {k + l)n. 

In this case, for all 1 < m < (A; + l)n, Vu{S) depends only on gm and /m-i in (3.19). Then 
the conclusion of Theorem 1.10 follows. 

4 Application examples 

In this section we give some application examples based on the examples given in [2]. 
Consider a Kukles type system of the form 

x = y, y = -{anxy + ao2|/^ + agox^ + aaix^y + ai2xy'^ + aosy^). (4.1) 

The authors [2J proved that if a^Q > and af^ — Sa^Q < then for (4.1) f 2 = f 4 = f e = 
f 8 = if and only if 021 = aos = aiiao2 = 0, which implies that the origin is a center. 
Moreover, there can be 3 limit cycles near the origin. See Theorem 4.1 in [2j and its proof. 
Based on this conclusion and by Theorem 1.6 we have immediately 
Proposition 4.1. Lei an, ao2, ctso, 021, ai2 and ao3 be bounded parameters satisfying 

aso > 0, - 8a3o < 0, |a2i| + jaosl + |aiiao2| > 0. 

Then there exists a neighborhood V of the origin such that the system (4.1) has at most 3 
limit cycles in V . 
Then consider 

X = -y + Ax"^ + Bxy + Cy^, y = x^ + xy'^ + y^. (4.2) 
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By Theorem 4.2 in [2] and its proof if < 2 then the origin of (4.2) is always a focus 
with \v2\ + Iv^l + \vq\ + \vs\ > 0. Moreover, there are systems inside (4.2) with at least 3 
limit cycles around the origin. Then by Theorem 1.6 again we have 

Proposition 4.2. Let A,B and C be bounded parameters with < 2. Then there 
exists a neighborhood V of the origin such that the system (4.2) has at most 3 limit cycles 
in V. 

Finally, consider 

k 

x = y, y = -{x^ + x^)-J2 b2jx^^y, (4.3) 

i=o 

where k > 2. By Theorems 1.7-1.9, we obtain 

Proposition 4.3. Let b2j be bounded parameters. Then 

(1) If bo = 0, the system (4.3) has at most k — 1 limit cycles near the origin; and 
k — 1 limit cycles can appear. 

(2) // 6o 7^ 0, there are systems inside (4.3) which have at least k limit cycles near 
the origin. 

References 

[1] Alvarez, M. J. & Gasull, A. [ 2005] Monodromy and stablility for nipotent critical 
points. Int. J. Bifurcation and Chaos 15, 1253-1265. 

[2] Alvarez, M. J. & Gasull, A. [ 2006] Generating limits cycles from a nilpotent critical 
point via normal forms, J. Math. Anal. Appl. 318, 271-287. 

[3] Amel'kin, V. V., Lukashivich, N. A. & Sadovski, A. P. [1982] Nonhnear Oscillations 
in Second Order Systems, Beloruss. Gos. Univ. Press, Minsk (in Russian). 

[4] Andreev, A. F. [1958] Investigation of the behavior of the integral curves of a system 
of two differential equations in the neighbourhood of a singular point, Transl. Amer. 
Math. Soc. 8, 183-207. 



25 



[5] Andreev, A. F., Sadovskii, A. P. & Tsikalyuk, V. A. [2003] The center-focus problem 
for a system with homogeneous nonhnearities in the case of zero eigenvalues of the 
hnear part, Diffenrential Equations 39, 155-164. 

[6] Bautin, N.N. [1952] On the number of hmit cycles which appear with the variation of 

coefficients from an equilibrium position of focus or center type. Mat. Sbornik N. S. 
30 181-196 (in Rassian); [1954] Translations Amer. Math. Soc. 100 (1954) 181-196 
(English translation). 

[7] Chavarriga, J., Giacomini, H., Gine, J. & Llibre, J. [2003] Local analytic integrabihty 
for nilpotent centers. Erg. Th. Dyn. Syst. 23, 417-428. 

[8] Chicone, C. and Jacobs, M. [1989] Bifurcation of critical periods for plane vector 
fields. Trans. Amer. Math. Soc. 312 433-486. 

[9] Han, M. [2000], The Hopf cycficity of Lienard systems, Appl. Math. Lett. 14(2), 
183-188. 

[10] Han, M. [2001], On some properties and limit cycles of Lienard systems. Dynamical 
systems and differential equations (Kennesaw, GA, 2000). Discrete Contin. Dynam. 
Systems, Added Volume, 426-434. 

[11] Liu, Y. & Li, J. [2009a] New study on the center problem and bifurcations of limit 
cycles for the Lyapunov system (I), Int. J. Bifurcation and Chaos 19, 3791-3801. 

[12] Liu, Y. & Li, J. [2009b] New study on the center problem and bifurcations of limit 
cycles for the Lyapunov system (II), Int. J. Bifurcation and Chaos 19, 3087-3099. 

[13] Liu, Y. & Li, J. [2010] Bifurcations of limit cycles and center problem for a class of 
cubic nilpotent system. Int. J. Bifurcation and Chaos 20(8), 2579-2584. 

[14] Liu, Y. & Li, J. [2011] Bifurcations of limit cycles created by a multiple nilpotent 
critical point of planar dynamical systems. Int. J. Bifurcation and Chaos 21(2), 1-8. 



26 



[15] Liapunov, A. M. [1983] Studies of one special case of the problem of stability of 
motion, Mat. Sb. 17, 253-333 (in Russian); [1966] Stability of Motion, Mathematics in 
Science and Engineering, Vol. 30, Academic Press, NY, London (Enghsh translation). 

[16] Moussu, R. [1982] Symetrie et forme normale des centres et foyers degeneres. Erg. 
Th. Dyn. Syst. 2, 241-251. 

[17] Romanovski V. G. & Shafer D. S. [2009] The Center and Cyclicity Problems: a 
Computational Algebra Approach. Birkhauser, Boston, MA. 

[18] R. Roussarie. [1998] Bifurcations of Planar Vector Fields and Hilbert's Sixteenth 
Problem. Progress in Mathematics, Vol. 164. Birkhauser, Basel. 

[19] Sadovskii, A. P. [1969] On the problem of distinguishing of a center and a focus. 
Differencial' nye Uravnenija. 5, 326-330 (in Russian). 

[20] Strozyna, E. & Zol^dek, H. [2002] The analytic normal for the nilpotent singularty, 
J. Diff. Eqs. 179, 479-537. 

[21] Takens, F. [1974] Singularties of vector fields, Inst. Hautes Etudes Sci. Publ. Math. 
43, 47-100. 

[22] Ye, Y. et al. [1986] Theory of Limit Cycles, Transl. Math. Monogr., vol.66, Amer. 
Math. Soc, Providence, RI. 



27 



